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Abstract 

A low momentum nucleon-nucleon (NN) potential Vi ow -k is derived from 
meson exhange potentials by integrating out the model dependent high mo- 
mentum modes of Vnn- The smooth and approximately unique Vi ow _k is 
used as input for shell model calculations instead of the usual Brueckner G 
matrix. Such an approach eliminates the nuclear mass dependence of the in- 
put interaction one finds in the G matrix approach, allowing the same input 
interaction to be used in different nuclear regions. Shell model calculations 
of 18 and 134 Te using the same input V/ otu _fc have been performed. For 
cut-off momentum A in the vicinity of 2 /m _1 , our calculated low-lying spec- 
tra for these nuclei are in good agreement with experiments, and are weakly 

dependent on A. 

21.60.Cs; 21.30.Fe; 27.80.+j 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



A fundamental problem in nuclear physics has been the determination of the effective 
nucleon-nucleon (NN) interaction used in the nuclear shell model, which has been success- 
ful in describing a variety of nuclear properties. There have been a number of successful 
approaches for this determination, ranging from empirical fits of experimental data, 

to deriving it microscopically from the bare NN potential. Despite impressive quantitative 
successes, the traditional microscopic approach suffers the fate of being "model dependent" 
owing to the fact that there is no unique Vnn to start from. Moreover, as the Brueckner G 
matrix has traditionally been the starting point, one obtains different input interactions for 
nuclei in different mass regions as a result of the Pauli blocking operator. 

In this Letter, we propose a different approach to shell model effective interactions that is 
motivated by the recent applications of effective field theory (EFT) and the renormalization 
group (RG) to low energy nuclear systems Our aim is to remove some of the model 

dependence that arises at short distances in the various Vnn models, and also to eliminate 
the mass dependence one finds in the G matrix approach, thus allowing the same interaction 
to be used in different nuclear regions such as 18 and 134 Te. A central theme of the RG- 
EFT approach is that physics in the infrared region is insensitive to the details of the short 
distance dynamics. One can therefore have infinitely many theories that differ substantially 
at small distances, but still give the same low energy physics if they possess the same 
symmetries and the "correct" long- wavelength structure |5|J§. The fact that the various 
meson models for Vnn share the same one pion tail, but differ significantly in how they 
treat the shorter distance pieces illustrates this explicitly as they give the same phase shifts 
and deuteron binding energy. In RG language, the short distance pieces of Vnn are like 
irrelevant operators since their detailed form can not be resolved from low energy data. 

Motivated by these observations, we would like to derive a low-momentum NN potential 
Vi ow -k by integrating out the high momentum components of different models of Vnn in 
the sense of the RG [§|§, and investigate its suitability of being used directly as a model 
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independent effective interaction for shell model calculations. We shall use in the present 
work the CD-Bonn NN potential |9j for Vjvjv- In the following, we shall first describe our 
method for carrying out the high- momentum integration. Shell model calculations for 18 O 
and 134 Te using V\ ow - k will then be performed. Our results will be discussed, especially 
about their dependence on the cut-off momentum A. 

The first step in our method is to integrate out the model dependent high momentum 
components of Vnn- In accordance with the general definition of a renormalization group 
transformation, the decimation must be such that low energy observables calculated in 
the full theory are exactly preserved by the effective theory. We turn to the model space 
methods of nuclear structure theory for guidance, as there has been much work in recent 
years discussing their similarity to the Wilson RG approach P JT0|JTT[| . While the technical 



details differ, both approaches attempt to thin-out, or limit the degrees of freedom one must 
explicitly consider to describe the physics in some low energy regime. Once the relevant low 
energy modes are identified, all remaining modes or states are "integrated" out. Their effects 
are then implicitly buried inside the effective interaction in a manner that leaves the low 
energy observables invariant. One successful model-space reduction method is the Kuo-Lee- 
Ratcliff (KLR) folded diagram theory |T2],[T3]]. For the nucleon-nucleon problem in vacuum, 



the RG approach simply means that the low momentum T matrix and the deuteron binding 
energy calculated from Vnn must be reproduced by Vi ow _ k , but with all loop integrals cut 
off at some A. Therefore, we start from the half-on-shell T-matrix 

T(k', k, k 2 ) = V NN (k', k) + j Q q 2 dqV NN (k>, q) fe2 _ g2 + . Q+ T(q, k, k 2 ). (1) 

We then define an effective low-momentum T-matrix by 

f A 1 

Ti ow - k (p',p,p 2 ) = V low -k(p',p) + / q 2 dqV low „ k (p',q)— T low _ k (q,p,p 2 ), (2) 

Jo p z — q z + t0 + 

where A denotes a momentum space cut-off (such as A=2/m _1 ) and {p',p) < A. We require 
the above T-matrices satisfying the condition 

T(p',p,p 2 ) = T low _ k (p\p,p 2 ); (p\p) < A. (3) 



The above equations define the effective low momentum interaction Vi ow -k- In the following, 
let us show that the above equations are satisfied by the solution 

V low . k = Q-Q' JQ + Q' J Q J Q-Q> J Q J Q J Q+ ... , (4) 

which is just the KLR folded-diagram effective interaction [r2JT5 |. A preliminary account 
of this result has been reported as a work in progress at a recent conference |14[ . 



In time dependent formulation, the T-matrix of Eq.(l) can be written as (k' 
VU(0, — oo) | k), U being the time evolution operator. In this way we can readily per- 
form a diagrammatic analysis of the T-matrix. A general term of it may be written as 

(*' I ( y + V W) V + V W) V W) v + • • •) I fc ) where < k ) = ( p - H °)i R o bein g the un P er " 

rurbed Hamiltonian. Note that the intermediate states (represented by 1 in the numerator) 
cover the entire space, and 1 = P + Q where P denotes the model space (momentum < A) 
and Q its complement. Expanding it out in terms of P and Q, a typical term of T is of the 
form V^-V^-V^V &V^V. Let us define a Q-box as Q = V + V^V + V^V^-V + ■■■, where 

eeeee 6 1 e ee 

all intermediate states belong to Q. One readily sees that the T-matrix can be regrouped as 
a Q-box series, namely (p' \ T \ p) — (p' | [Q + Q^Q + Q^Q^Q + • • •] | p). Note that all 
the Q-boxes have the same energy variable, namely p 2 . 

This regrouping is depicted in Fig. 1, where each Q-box is denoted by a circle and the 
solid line represents the propagator — . The diagrams A, B and C are respectively the one- 
and two- and three-Q-box terms of T, and clearly T=A+B+C+- • •. Note the dashed vertical 
line is not a propagator; it is just a "ghost" line to indicate the external indices. We now 
perform a folded-diagram factorization for the T-matrix, following closely the KLR folded- 



diagram method |p!z[^3[ . Diagram B of Fig. 1 is factorized into the product of two parts 
(see Bl) where the time integrations of the two parts are independent from each other, each 
integrating from — oo to 0. In this way we have introduced a time-incorrect contribution 
which must be corrected. In other words B is not equal to Bl, rather it is equal to Bl 
plus the folded-diagram correction B2. Note that the integral sign represents a generalized 



folding [12,13 



Similarly we factorize the three-Q-box term C as shown in the third line of Fig. 1. 
Higher-order Q-box terms are also factorized following the same folded-diagram procedure. 
Let us now collecting terms in the figure in a "slanted" way. The sum of terms Al, B2, C3... 
is just the low-momentum effective interaction of Eq.(4). (Note that the leading Q-box of 
any folded term must be at least second order in Vat at, and hence it is denoted as Q'-box 
which equals to Q-box with terms first-order in Vnn subtracted.) The sum Bl, C2, D3.... 
is VJ oa ,_fc— Q. Similarly the sum C1+D2+E3+- ■ ■ is just Vi ow _k— Q^Q- (Note diagrams Dl, 
D2, ■ • -, El, E2, • ■ ■ are not shown in the figure.) Continuing this way, it is easy to see that 
Eqs. (1) to (3) are satisfied by the low momentum effective interaction of Eq.(4). 

The effective interaction of Eq.(4) can be calculated using iteration methods. A number 
of such iteration methods have been developed; the Krenciglowa-Kuo [15| and the Lee-Suzuki 



iteration methods |16| are two examples. These methods were formulated primarily for the 
case of degenerate PH P, H being the unperturbed Hamiltonian. For our present two- 
nucleon problem, PHqP is obviously non-degenerate. Non-degenerate iteration methods 



17| are more complicated. However, a recent iteration method developed by Andreozzi 18 



is particularly efficient for the non-degenerate case. This method shall be referred to as the 
Andreozzi-Lee-Suzuki (ALS) iteration method, and has been employed in the present work. 

We have carried out numerical checks to ensure that certain low-energy physics of Vnn 
are indeed preserved by Vi ow _k- We first check the deuteron binding energy BE^ given by 
Vi ow -k. For a range of A, such as 0.5/m -1 < A < 3/m _1 , BEd given by Vi ow -k agrees very 
accurately (to 4 places after the decimal) with that given by Vnn- in Fig. 2, we present 
some 1 5'o phase shifts calculated from the CD-Bonn Vnn (dotted line) and the Vi ow _k (circles) 
derived from it, using a momentum cut-off A = 2.0/m -1 . As seen, the phase shifts from the 
former are well reproduced by the latter. We have also checked the half-on-shell T-matrix 



given by Vnn and by V ow -k, and found very good agreement between them fL4| . In short, 
our numerical checks have reaffirmed that the deuteron binding energy, low energy phase 
shifts and low momentum half-on-shell T-matrix of Vnn are all preserved by V ow -k- As far 
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as those physical quantities are concerned, Vi ow -k and Vnn are equivalent. 

Having proven the "physical equivalence" of Vi ow -k and Vnn in the sense of the RG, 
we turn now to microscopic shell model calculations in which we use Vi ow -k as the input 
interaction. We have performed shell-model calculations for 18 and 134 Te following the same 
procedure as outlined in Refs. , except that the G-matrix vertices used there are replaced 
by our present Vi ow -k- A model space with two valence nucleons in the (0^5/2,0^3/2, lsi/2) 
shell is used for 18 0, and a (0^7/2, 1^5/2, 1^3/2, 2si/ 2 , 0h u / 2 ) one for 134 Te. In a concurrent 
paper fll]] we have found that Vi ow -k is almost independent of the underlying Vnn for the 



values of A considered here. Therefore, although the CD-Bonn potential [[J is used in our 
calculations, we stress that very similar results will be obtained if we calculate Vi ow -k from 
other models such as the Paris or Argonne V-18 potentials. 

Our calculated low-lying J n states of 18 and 134 Te are presented in Fig.3 and Fig.4. 
In the same figures, results of the corresponding G-matrix calculations are also shown; 
the Vi ow _k results are just as good or slightly better. It may be mentioned that the G- 
matrix is energy dependent and Pauli blocking dependent, while V\ ow -k is not which is a 
desirable feature. An important issue is what value one should use for A. Guided by general 
EFT arguments, the minimum value for A must be large enough so that V ow ^k explicitly 
contains the necessary degrees of freedom for the physical system. Such a value for A min is 
signalled when the calculated spectra first become insensitive to A H . Conversely, we want 
A to be smaller than the short distance scale A max at which the model dependence of the 
different Vnn starts to creep in ||. Systems in which these two constraints are consistent 
with each other (i.e., A min < A max ) are amenable to EFT-RG inspired effective theories, as 
they possess a clear separation of scales between the relevant long wavelength modes and the 



model dependent short distance structure. We have found JTT|] that A max should not be much 
greater than 2.0-2.5 fm _1 as this is the scale at which V ow _k first becomes dependent on the 
particular Vnn used. There is another consideration: Most NN potentials are constructed 
to fit empirical phase shifts up to Ei ab « 350 MeV 0. Since Ei ab < 2h 2 A 2 /M, M being the 
nucleon mass, and one should require V ow -k to reproduce the same empirical phase shifts, 
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a choice of A in the vicinity of 2 fm~ l would seem to be appropriate . 

Guided by the above considerations, we have used in our calculation two values for the 
momentum cut-off, namely A = 2.0 and 2.2 fm _1 . It is satisfying to see that the results in 
Fig. 3 and Fig. 4 are rather insensitive to the choice of A; this suggests that A rnin is around 
2fm~ 1 in harmony with the EFT philosophy mentioned earlier. Perhaps more importantly, 
both are in satisfactory agreement with experiments. Note that all energies are relative to 
the ground state energy of the respective core nucleus. 

We emphasize that we have used the same Vi ow -k interaction in both 18 and 134 Te 
calculations, and it appears to work equally well for both nuclei. This is in marked contrast 
to the traditional approach in which one uses G-matrices, as the Pauli blocking operator 
is very different for the two nuclear mass regions. This is an exciting result, as it suggests 
the possibility for a common shell-model interaction that is nearly model independent and 
suitable for a wide range of nuclei. 

In summary, we have investigated a RG-EFT inspired approach to shell model calcula- 
tions that is a "first step" towards a model independent calculation that uses one common 
interaction over a wide range of nuclei. Using the KLR folded diagram approach in conjunc- 
tion with the ALS iteration method, we have performed a RG decimation where the model 
dependent pieces of Vnn models are integrated out to obtain a nearly unique low momentum 
potential V[ ow _k- This Vi ow _k preserves the deuteron pole as well as the low energy phase 
shifts and half-on-shell T matrix. We have used Vi ow ^ k , which is a smooth potential, directly 
in shell model calculations of 18 and 134 Te without first calculating the G matrix. The re- 
sults are in satisfactory agreement with experiment for both nuclei, and they are insensitive 
to A in the neighborhood of A 2 fmr x . For completeness, we add that satisfactory results 
have also been obtained using the same V[ ow -k in shell model calculations in the 208 Pb and 
132 Sn regions. Calculations of spin orbit splittings using the same V/ ow _fc have also led to 
satifactory results. These results are a work in progress and will be reported in a subsequent 
paper. We do feel that Vi ow -h may become a promising and reliable effective interaction for 
shell model calculations of few valence nucleons, over a wide range of nuclear regions. 
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FIG. 1. Folded-diagram factorization of the half-on-shell T-matrix. 
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FIG. 2. Comparison of phase shifts given by V/ ow _fc and Vnn- 



11 



-8.0 _ 



-9.0 _ 



> 

a; 



-10.0 _ 



-11.0 _ 



-12.0 _ 





18q 






2+ 


2+ 


2+ 


2+ 
0+ 


o+ 


n+ - 


0+ 




4+ 


4 + 


4+ 


4+ 


2+ 




2+ 


2+ 


— 


2+ - 




0+ 


0+ 


0+ 


0+ 




Expt. 


A = 2.0 A 


= 2.2 


G matrix 



FIG. 3. Low lying states of 18 O. 
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FIG. 4. Low lying states of 1M Te. 
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